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We discuss the possibility of a quantum phase transition in ultracold spin polarized fermionic gases 
which exhibit a p-wave Feshbach resonace. We show that when fermionic atoms form a condensate 
that can be externally tuned between the BCS and BEC limits, the zero temperature compressibility 
and the spin susceptibility of the fermionic gas are non-analytic functions of the two-body bound 
state energy. This non-analyticity is due to a massive rearrangement of the momentum distribution 
in the ground state of the system. Furthermore, we show that the low temperature superfluid density 
is also non-analytic, and exibits a dramatic change in behavior when the critical value of the bound 
state energy is crossed. 
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Introduction: Recent experiments in cold fermionic 
gases have shown that s-wave magnetic field induced Fes- 
hbach resonances can be used to form diatomic molecules 
of 40 K 1] and 6 Li 0,H||, 

which undergoe Bose-Einstein 
condensation (BEC) on the higher magnetic field side of 
the resonance. On the lower magnetic field side of the 
resonance, it has also been established that Cooper pair- 
ing takes place and a BCS condensate is formed. These 
studies in cold fermionic gases led to the first experi- 
mental realization of the theoretically proposed BCS-to- 
BEC crossover in three dimensional continuum s-wave 
superfluids 0, Three early theoretical works that 
considered the possibility of s-wave superfluidity in the 
context of (what is known today as) the BCS-to-BEC 
crossover should be highlighted. The first is by Ea- 
gles 0, where the possibility of pairing without con- 
densation is described in a continuum model in the con- 
text of superconductors with low carrier concentration. 
The second is Leggett's seminal work 8], in which the 
T = s and p-wave BCS-to-BEC evolution are dis- 
cussed as a crossover phenomenon in the context of a 
variational ground state wavefunction. And the third is 
the work of Nozieres-Schmitt-Rink where the s-wave 
BCS-to-BEC crossover in a lattice is described. Further- 
more, much of the theoretical [lol ITl| and experimental 
01111 efforts that followed described only the BCS- 
to-BEC crossover in s-wave systems. 

In this manuscript, we present a functional integral 
analysis of the BCS-to-BEC evolution in p-wave fully 
spin-polarized Fermi gases, where p-wave Feshbach res- 
onances have already been observed 0, ^| . We show 
that a quantum phase transition takes place when the 
chemical potential crosses a critical value, instead of the 
usual smooth BCS-to-BEC crossover that occurs in s- 
wave superfluids [l4j |. The atomic compressibility and 
the spin-susceptibility of the Fermi gas are computed 
and are shown to be non-analytic in the p-wave case, 
as a consequence of a major rearrangement in the mo- 
mentum distribution as the critical point is approached. 
This non-analytic behavior suggests the occurrence of a 
quantum phase transition, which is further confirmed by 
a discontinuous change in the temperature dependence of 



the superfluid density of the gas at the transition point, 
which goes from power-law on the BCS side of the reso- 
nance to exponential on the BEC side of the resonance. 

We study the case of two-dimensional systems, which 
can be prepared experimentally through the formation 
of a one-dimensional optical lattice, where tunnelling be- 
tween lattice sites is suppressed by a large trapping po- 
tential. The form of the trapping potential can be chosen 
to be Vtrap = — /o exp[— 2(x 2 + y 2 )/w 2 ] cos 2 (fc z z), where 
2n/k z is the wavelength of the light used in the laser 
beam. We assume that the width w is such that w ^> \p, 
where Xf = 27r/fcp is proportional to the interparticle 
spacing of a Fermi gas with Fermi wavevector k^?, such 
that the problem is essentially two-dimensional. 

Hamiltonian: We study a two-dimensional continuum 
model of spin polarized (all atoms in the same hyperfine 
state) fermionic atoms of mass m and density n — kp /4ir. 
In the presence of an external magnetic field h, the sys- 
tem is described by the Hamiltonian (H = ks = 1) 

W = E^kT^T + \ E ^k4<A'q: (1) 
k k,k',q 

where fr kq = V , -k+q/2TV'k+q/2T and 6c = £k - P, with 
ek = k 2 /2m, and p, = fi + ga^B^i- The direction of 
the magnetic field h, which was chosen to define the spin 
quantization axis z, need not to coincide with the spatial 
direction z of propagation of the laser beam. 

The interaction potential is approximated by the fol- 
lowing separable function in k-space, 

Vkk' = -Ar(k)r(k'), (2) 

where A is the interaction strength and T(k) = 
h(k)cos((p), where the funcion h(k) = (k/ki)/[l + 
k/ko] 3 / 2 controls the range of the interaction, <p is the 
momentum angle in polar coordinates, and Rq ~ fc^ 1 
plays the role of the interaction range. The functional 
form of r(k) can be shown to produce the correct asymp- 
totic behavior at small and large momenta |l4j . and 
its angular dependence reflects equal contributions from 
the angular momentum channels I = ±1. In the limit 



2 



of small momenta, this approach is identical to the T- 
matrix formalism but has the added advantage of 
making unnecessary to introduce the scattering length 
as a relevant parameter, which is quite problematic in 
two-dimensions The BCS-BEC evolution can be 

safely analyzed provided that the system is dilute enough 
{kp -C fcg), i.e., the square of the interparticle spacing 
(~ fc^ 1 ) is much larger than the square of the interaction 
range (~ k^ 1 ). Throughout the manuscript, we choose 
to scale all energies with respect to the Fermi energy 
£f = kp/2m and all momenta with respect to k f- 

Effective Action: The partition function Z at a temper- 
ature T — is written as an imaginary-time functional 

integral with action S = dr[J2^ V'kt CO^W'kTC'") +^]- 
Introducing the usual Hubbard-Stratonovich field 0q(V), 
which couples to ^T^t ; anc i integrating out the fermionic 
degrees of freedom, we obtain 



Z = J V(t>V(j)* exp(-S* cff [ 



(3) 



with the effective action given by 



S e ff = 



(It 



U(t) 



E 

k,k' 



£k j- 

2" k,k ' 



Trln 2 G k,k' ( r ) 



where U(t) = £ k |^ k (r)| 2 /(2A) and G kk ,(r) is the (in- 
verse) Nambu matrix, 



Gt.Ur) = 



~{d T +Ck)^k,k' A k>k /(r) 



A^, k (r) 



-(dr ~ Ck)^k,k' 



(4) 



with A k)k ,(r) - k _ k ,(r)r((k + k')/2). 

Saddle Point Equation: After Fourier transforming 
from imaginary time to Matsubara frequency {ik n — 
i(2n + 1)tt//3) and performing the frequency sum, the 
saddle point condition [SSeff/S^)* (t')]a = can be cast 
in the form of the familiar order parameter equation, 



i = E 



^tanhf^ 
2E k V 2 



(5) 
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FIG. 1: Universal plot (for any magnetic field hg) of jl = 
fi + gufiBhz and order parameter amplitude Ao (inset) as 
functions of Eb = Eb(hz) + 2g^UBhz for ko = fci = 10kp in 
the spin polarized p-wave case. All quantities are in units of 



number equation as Ao = X) k n k, where the momentum 
distribution n k is given by 



1 

n k = - 
2 



1 — tanh f — — - 
Ev \ 2 



(7) 



Thus, at T = 0, the saddle point and number equations 
reduce to E k r2 ( k )[( 2e k - Eb^ 1 - (2£k) -1 ] = and 
N = X)k(l — Ck/^k)/2, respectively. The solutions for 
Ao and p, at T = as functions of the binding energy 
Et in the case of p-wave pairing symmetry are plotted in 
Fig. 1 for fco = ki = 10k p. The point jj, = is achieved 
for Eb = — 1.087eF and corresponds to A = 19.063e_F. 

Gaussian Fluctuations: We now investigate the effect 
of Gaussian fluctuations in the pairing field 4>q (r) about 
the static saddle point value Ao- Assuming (j)q(T) = 
Ao^q.o + ?7q(T) and performing an expansion in S e s to 
quadratic order in 77, one obtains 



S Gauss = S [A ] + ^ n }(q)M(q)ri(q), 



(8) 



where i? k 



+ |A k | 2 is the quasiparticle excitation 



energy, and A k = AoT(k) plays the role of the order pa- 
rameter function. We eliminate the interaction strength 
A in favor of the two-body bound state energy Eb(hz) 
in vacuum (and in the presence of a magnetic field) by 
using the relation 1/A = r 2 (k)/(2e k — Eb), where 
Eb = Eb(hz) + 2gzi/J,Bh z . The renormalized gap equa- 
tion in terms of Eb then takes the form 



1 



tanh(/3£; k /2) 



2e k - E b 



2Eu 



= 0. (6) 



Number Equation: Using the relation N = —dQ/dj-i 
and the saddle point approximation for the thermody- 
namic potential, £Iq — S c s [Aq]//9, one can write the 



where Sq is the saddle point action, the vector n(q) is 
such that rf(q) — [i]*(q),i](—q)], and q = (q,iq m ), where 
iqm = i2mn//3 is a bosonic Matsubara frequency. The 
2x2 matrix M(g) is the inverse fluctuation propagator. 

The Gaussian fluctuation term in the effective action 
leads to a correction to the thermodynamic potential, 
which can be rewritten as f^Gauss = fio + ^fluct, with 
f^fluct = P^ 1 J2 q mdet[M(<7)]. Therefore, using the rela- 
tion N = —dQ/dfj,, one can write the corrected number 
equation as Na anss = N + Nf[ uct , where N is the saddle- 
point level number of particles given above, and 



AWt = 



on 



fluct 



dfi 



T EE 



-d(dctM)/d M 
detM(q, iq n ) 



(9) 
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At low T, the Goldstone mode u> — c|q| dominates the 
contribution to iVfl uct , leading to 

T 2 A r)r 

IVfluct -"^C(3)4^ 3 , (10) 

27T C™ 6[l 

which vanishes in the limit of T — > 0. Therefore, analo- 
gously to the three-dimensional s-wave case @, Eq. (7) 
provides a very accurate description of the number equa- 
tion near and at T = 0, thus confirming Leggett's sug- 
gestion However, it is well known that the same is 
not true near T c , where the effects of temporal fluctua- 
tions are essential to describe the BEC regime 5j. The 
discussion of this interesting limit will be postponed to 
a future manuscript, and we will focus here on the low 
temperature properties, to be discussed next. 

Momentum distribution: The momentum distribution 
rik given by Eq. (7), which at zero temperature reduces 
to nk = (1 — £k/-Ek)/2, i s plotted in Fig. 2 for the case 
of p- wave pairing symmetry as a function of k = (k x , k y ), 
together with the contour plots. Notice that becomes 
discontinuous when the chemical potential crosses zero, 
which coincides with the collapse of the two Dirac points 
to a single point k = and the appearence of a full 
gap in the quasiparticle excitation spectrum. This ma- 
jor rearrangement of the momentum distribution has a 
dramatic effect in the atomic compressibility, which is 
discussed next. 

Atomic Compressibility: The first derivative of the 
chemical potential with respect to the density n = N/ L 2 
becomes non-analytic at the critical value of the bind- 
ing energy in the p-wave case. As a consequence, the 
isothermal atomic compressibility k, defined by 

_ L 2 d 2 n _ 1 dn 

will develop a cusp when expressed in terms of Eb, its 
first derivative with respect to Eb diverging at the criti- 
cal point, as shown in Fig. 3. In the s-wave case, however, 
k is smooth for all values of Eb 14] . This non-analytic 
behavior of the p-wave atomic compressibility, combined 
with the appearence of a full gap in the excitation spec- 
trum, suggests the existence of a quantum critical point 
at jl = 0. 
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FIG. 3: Plot of dn/d^, (in units of k%/AirtF) and its first 
derivative with respect to Ef, (inset) as functions of Eb in the 
case of spin polarized p-wave pairing and ko — ki = IO/cf. 



Spin Susceptibility: The phase transition discussed in 
the previous section also manifests itself in the spin sus- 
ceptibility. The application of a small probe magnetic 
field Hz along the same direction (z) of h generates the 
spin susceptibility response Xzz = (—l/L 2 )(d 2 fl/dH 2 ), 
which can be rewritten in the case of spin polarized atoms 
as 

lo 2 9 2 ^l 2 2 dn 
X^ = - T2 9^b^=9^b-^- (12) 

Thus, the graph in Fig. 3 also represents a universal plot 
°f Xzz I ' glik^B as a function of Eb- 

Superfluid Density: We now turn our attention to the 
behavior of the low temperature superfluid dentity ten- 
sor pij(T,Eb) as the critical value of the binding en- 
ergy Eb is crossed. This tensor is associated with phase 
twists of the superconductor order parameter [T^ . and 
can be obtained by taking </> q —> </> q exp(z(? q ) and ex- 
panding the effective action S e g in powers of q about 
the saddle point with # q = 0. The resulting difference 
in the action, AS = S' c ft(0 q ) — S e s(9q = 0), becomes 
AS{T) = -{L 2 /2)E q V-q^S^ijCr), with the super- 



4 



-0.005 - 



-0.01 



-0.015 - 



-0.02 



-0.025 






1 1 1 1 l_ 




" 



0.0005 
T 



0.001 



FIG. 4: Plot of Ap(T)/T 2 (in units of e~ 2 ) as a function of 
temperature (in units of sf) for various values of the binding 
energy E},. Inset: Zero-temperature slope of Ap(T)/T 2 (in 
units of eZ 3 ) as a funcion of Eh. 



fluid density tensor given by 

Pij (T) = ^ E Pnk^fliCk - Ykdi&dj^] , (13) 

k 

where 7ik is the momentum distribution, Ik = 
(2T)" 1 sech 2 (£; k /2r) is the Yoshida distribution, and 
denotes the partial derivative with respect to fcj. Notice 
that = p yy = p, while p xy — p yx = 0. In addi- 
tion, notice that at T — 0, Pij(0) = n/m, such that 
dpij/dp = (l/m)dn/dfx and dpij/dH^ = {l/m)dn/dHz. 
Using our energy and momentum scales, we define the 
dimensionless quantity Ap(T) = mp(T)/n — 1, which is 
shown in Fig. 4 as a function of temperature for differ- 
ent values of the binding energy. The linear behavior 
of Ap{T)/T 2 for values of Eb that correspond to ft > 
indicates a T 3 dependence of the superfluid density on 
temperature on the BCS side of the transition. This be- 
havior is in fact confirmed by our analytical calculation 



of Ap(T) at low temperatures and in the case of short 
range interactions (fco — * oo). In the BCS limit, we found 
Ap(T) ~ CT 3 , with the coefficient C weakly dependent 
on Eb- This power-law behavior reflects the nodal (gap- 
less) structure of the p-wave excitation spectrum. In the 
BEC limit, we obtained Ap(T) ~ exp(-|/2|/T), the ex- 
ponential behavior reflecting the appearance of a full gap 
to the addition of quasiparticles for p, < 0. Fig. 4 also 
shows (inset) the zero temperature slope of Ap(T)/T 2 as 
a function of the binding energy Eb, which is clearly dis- 
continuous at the critical point Eb = — 1.087ej?. These 
results further confirm the existence of a quantum phase 
transition along the BCS-to-BEC evolution as a function 
of interaction strength (binding energy) in the case of 
p-wave spin polarized atoms. 

Summary: We proposed the existence of a quantum 
phase transition in the BCS-to-BEC evolution of p-w&ve 
fully spin polarized Fermi gases as a function of the two- 
body bound state energy. We have shown that, at a 
critical value of this binding energy, the momentum dis- 
tribution undergoes a major rearrangement in k-space, 
which leads to a non-analytic behavior of the atomic 
compressibility and spin susceptibility of the gas. Fur- 
thermore, the low temperature superfluid density of the 
system presents a dramatic change in behavior as the 
critical point is crossed, its temperature dependence go- 
ing discontinuously from power-law on the BCS side of 
the transition to exponential on the BEC side of the tran- 
sition. 

We conclude by suggesting that this phase transition 
may be observable in traps of 6 Li and 40 K g ases which 
exhibit p-wave Feshbach resonances |12| 13]. The oc- 
currence of this phase transition may be investigated 
through the direct measurement of the atomic compress- 
ibility, spin susceptibility or superfluid density as func- 
tions of binding energy or magnetic field. 
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